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SUMMARY 


Using anisotropic elasticity theory, delamination stresses in a semicircular laminated composite curved bar 
subjected to end forces and end moments were calculated, and their radial locations were determined. A 
family of design curves was presented, showing variation of the intensity of delamination stresses and their 
radial locations with different geometry and different degrees of anisotropy of the curved bar. The effect 
of anisotropy on the location of peak delamination stress was found to be small. 

INTRODUCTION 

The major cause of degradation in stiffness and strength of laminated composite materials is the growth 
of delamination between individual composite layers. Excess delamination may result in ultimate fatigue 
failure. The most common failure mode of laminated composite material is delamination failure. Delam- 
ination may result from the following: eccentricities in the structural load paths, inducing out-of- plane 
loads; discontinuities in the structure, creating local out-of-plane loads; low-velocity impacts; cyclic load- 
ing; incomplete curing; and the introduction of foreign particles during the manufacturing process. The 
delamination growth may redistribute stresses in the composite plies and therefore may reduce the residual 
stiffness and strength of the laminated composites. In the common application of the composite materials, 
many structural parts have certain degrees of curvature. Because of this curvature, radial stress can be 
generated. If the radial stress is in tension, this stress will function as open-mode delamination stress and 
will cause the delamination to initiate and grow under service loading (mostly cyclic loading). 

Delamination problems in composite materials have been studied extensively recently, and major re- 
search activities are listed in O’Brien (1984a). There are a variety of ways to conduct composite delam- 
ination studies using different types of test specimens. However, one of the most attractive test coupon 
geometries is in the shape of a semicircular curved bar. It is well known that if the curved bar is subjected 
to bending, radial stress and shear stress can be induced inside the curved bar. If the loading is intense 
enough, open-mode delamination can take place at the site of peak tensile radial stress. If the composite 
curved bar has weak shear strength, shear-mode delamination can occur at the site of peak shear stress. 
Since the peak value of tensile radial stress (or peak shear stress) occurs only at a particular point, the 
curved bar offers an excellent situation for studying the initiation and subsequent growth of delammation 
under cyclic loading and the fatigue behavior (degradation of stiffness and strength) of laminated composite 
materials. 

The purpose of this report is to document the calculation of delamination stress (open-mode or shear 
mode) and its exact radial location and to relate how delamination stress and its location change with the 
degree of anisotropy (different stacking sequences) and the wall thickness of the curved bar. 

NOMENCLATURE 

a inner radius of curved bar, in 

b outer radius of curved bar, in 

El modulus of elasticity of single ply in fiber direction, lb/in 2 

E t modulus of elasticity of laminated composite in r direction, lb/in 2 

£ t modulus of elasticity of single ply in direction transverse to fiber direction, lb/in 2 

Eg modulus of elasticity of laminated composite in 0 direction, lb/in 2 

e loading axis offset, in 

G\si shear modulus of single ply, lb/in 2 

Gg r shear modulus of laminated composite, lb/in 2 

h width of curved bar, in 


k anisotropic parameter, \JEgj E T 

M applied end moment, in-lb 

P applied end force, lb 

r radial coordinate, in 

r m radial location of (< 7 r ) max , in 

radial location of (a') msx , in 

t thickness of curved bar, t = (6 — a), in 

/? anisotropic parameter, f3 = i/l -f (Eg/E r )(l + 2i/ 6r ) -(- Eg/Ge r 

9 tangential coordinate, rad 

l 'LT>* / ‘rL Poisson ratios of single-layer composite 
vg r Poisson ratio of laminated composite 

&D maximum radial stress (open-mode delamination stress) induced by both 

P and Af, op (^v) m&x “f* (o r )max) lb/in 2 

<r r radial stress induced by end forces P, lb/in 2 

o' r radial stress induced by end moments M , lb/in 2 

(^r)max maximum radial stress induced by end forces P, lb/in 2 

(o’r)nuut maximum radial stress induced by end moments M, lb/in 2 

c 7 & tangential stress induced by end forces P, lb/in 2 

erf, tangential stress induced by end moments M, lb/in 2 

T r g shear stress induced by end forces P, lb/in 2 

T ' r6 shear stress induced by end moments M , lb/in 2 

COMPOSITE CURVED BAR 

Figure 1 shows a semicircular composite curved bar subjected to end forces P. Because finite area is 
needed for load attachment in the fatigue tests, both ends of the curved bar have to be extended slightly. 
Thus, the loading axis will have slight offset e from the vertical diameter of the curved bar. Therefore, the 
loading state in figure 1 (or fig. 2, A) is the combination of two cases: 

1. Bending due to end forces P with the loading axis coinciding with the vertical diameter of the curved 
bar (fig. 2, B) 

2. Pure bending due to end moments M(= Pe) created by the loading axis offset e (fig. 2, C) 

It has been observed that the highest probability of delamination onset takes place at the interfaces of 0° 
and 90° plies because of high Poisson’s ratio mismatch (O’Brien, 1982; O’Brien, 1984b). In constructing 
the curved-bar test coupon, it is desirable to introduce 90° plies at, or in the vicinity of, a peak stress 
point (that is, a peak radial tensile stress point or peak shear stress point) and thereby ensure that the 
delamination will initiate at the peak stress point, which is yet to be determined. In the following sections, 
the peak radial tensile stress (or peak shear stress) and its radial location will be calculated. 

ANISOTROPIC CURVED BAR UNDER END FORCES 

Figure 3 shows the anisotropic semicircular curved bar subjected to end forces P with the loading axis 
coinciding with the vertical diameter of the curved bar. The situation in figure 3 is similar to that in 
figure 2, B. 
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Stresses Induced 


If the composite material of the curved bar is treated as a continuous anisotropic material, then the 
stresses induced in the composite curved bar due to the end forces P may be written as (Lekhnitshii and 


others, 1968, p. 99, taking 

u = 0°) 

<Tr(r,6) = 

P b 

(r\& ( a\P /b\P 

bhgi r 

+ U) (?) 

oe{r\8) = 

P b 
bhgi r 

d+w^' + d- 

T r e{ri9) = 

P b 

/ r \P / a\ @ /6\ ^ 

bhgi r 

\ftj \b) \r/ 


.0 /n\P 


sin 8 




cos 8 


( 1 ) 

( 2 ) 

( 3 ) 


where 

a is inner radius of the curved bar, 
b outer radius of the curved bar, 

h width of the curved bar, 

r radial coordinate, 

$ tangential coordinate, 

<r r radial stress, 

a# tangential stress, 

T r $ shear stress, 

and 



and the anisotropic parameter (3 is defined as 


0 — 


' 1 + "—-(I ” 2 l/0 r ) + 

hj r 


Ee_ 
G e r 


( 4 ) 

( 5 ) 


where 

E$ is modulus of elasticity in 6 direction, 

E r modulus of elasticity in r direction, 

G$r shear modulus, and 
u$ r Poisson’s ratio. 

For isotropic materials, (3 = 2. 

Notice from equations (1) and (3) that the magnitudes of <7 r and T r $ are identical, but they are out of 
phase by 7r/2. The maximum value of <J r occurs at cross section 6 = x/2 (midspan of the curved bar), and 
T r e reaches its peak value at the two load application cross sections 0 = 0 and 6 - 7r: 

0 Tr (r, |) = -T re (r, 0) = T r9 (r, x) (6) 
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Thus, the semicircular curved-bar test coupon can provide the same intensities of open-mode and shear- 
mode delamination stresses simultaneously. If the composite is weak in open-mode strength, the delami- 
nation will initiate at the midspan 0 = 7t/ 2. On the other hand, if the composite is weak in shear strength, 
delamination will start at both ends of the curved bar (that is, 0 = 0 and 0 = ir). 

Because of the relationship between a r and r r g given in equation (6), analysis will be limited to a r . 
Equation (1) (for 0 = 7t/ 2) may be rewritten as 



P a 
haB r 


am*- cr- 


( 7 ) 


where 



Or, in dimensionless form, 


h(b- a) 

P 


<7 r 



( b/a — 1) a 


B 





(8) 


( 9 ) 


For the isotropic case (that is, /? — 2), equation (7) may be written in the following form (for h = 1): 


( 7T V 

P 

a 2 b 2 I., 9, 

° T \ 2 /. 

0=2~~N 

r + 3 (a + b 2 ) 

r 


( 10 ) 


where 

N = -a\B) 0=2 = a 2 - b 2 + (a 2 + & 2 )ln - (11) 

Notice that equation (10) is identical with the expression given by Timoshenko and Goodier (1970) for 
isotropic materials). 


Location of Maximum a r 


The radial location r = r m where a r (r, tt/ 2) reaches its peak value (<r r ) max may be found by differentiating 
equation (7) with respect to r and setting the resulting derivative to zero, or 


d_ 

dr 


0> 



= 0 


( 12 ) 


from which r m is found as 



1 

2(/*-l) 


■ 

fb \ 0 i 

2 fbV 

it 

|(i) +i J 

+ 4 (^ -,)(-) 



(13) 


For the isotropic case 




( 14 ) 
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((7 r )max may be written as 


h(b — a) 


(^r)max — 


h(b - a) 


Or T 


■•!) = 


( b/a — 1) a 
B r m 


(tulY + f!Wiy_f-)-i 

V a J \aj \r m J \aj 


(15) 


where r m is given by equation (13). 


Thin Curved Bar 


When the thickness t (= 6 - a) of the curved bar is small, 


- = 1 + - 
a a 



(16) 


The expression for B given in equation (8) may be expanded in terms of small quantity t/a as 


( B)i 




+ 


(17) 


Notice that the lower order terms up to Oit/a) 2 cancelled out. 

The expression for r m given in equation (13) may be expanded to the following form: 

2 


^ =1 + 

® / f> — 


1 i /tv 

2 a 4 \aj 


+ 


(18) 


Notice that r m is independent of (3 and that (cr r ) raax is always located on the inner side of the middle 
surface of the curved bar. 

Lastly, equation (15), with equation (18) applied, may be expanded as 


'h(b — a) 


Mr 


J 6- 


- _I (l\ 3 ?L + 

►a B\a) 4 + 


(19) 


where the lower order terms up to 0(t/a ) 2 cancelled out. 
Combining equations (17) and (19) there results 


h(b-a) 

P 



3 

2 


( 20 ) 


which is independent of anisotropy (that is, the value of /?). 

ANISOTROPIC CURVED BAR UNDER END MOMENTS 

Figure 4 shows the anisotropic semicircular curved bar under pure bending due to end moments M. This 
case corresponds to that in figure 2, C. 
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Stresses Induced 


If the composite material of the curved bar is treated as continuous anisotropic material, the radial stress 
<r r , tangential stress <7^, and shear stress t' t q induced in the curved bar under the end moments M may be 
expressed as (Lekhnitskii and others, 1968, p. 97) 


a' r {r) 


M 

1 _ 

1 - (a/b) k+1 / 

r\ k - 1 1 — (a/6)* _1 / 1 

a\ k+1 ( 

b\ k+1 ' 

b 2 hg 


1 — (a/b) 2k V 

6/ 1 — (a/6) 2 * V 

b) \ 

r) , 


°»( r ) 


M 

b 2 hg 


<6 = 0 


i iwr, 

1 - ( a/b) 2k 



l - fr / d)*- 1 / qN*- 1 
l-(a/b) 2k \b ) 


CT 


( 21 ) 


( 22 ) 

(23) 


where the anisotropic parameter k is defined by 


Wf 


and 


__ 1 - (a/*>) 2 k_ [l ~ (a/f>) fc+1 ] 2 i k(a/b) 2 [l - (a/6)*- 1 ] 

2 k+l [1- ( a/6 ) 2 * ] + k - 1 [1 - (a/6) 2*] 


Rewriting a' T {r ) (eq. (21)) in terms of 6/a instead of a/6 yields 


K(r) = ~ 


M 

a 2 hC 


![©”-] 


aj 


*+i 


- 1 


(r- fern ere 


-(*+!)' 


where 


HicHieH-rMcrn'-Acnt 


6 \ * 

,a) ^ 


In dimensionless form, equation (26) may be written as 


ha m (b - a) , 

— M — = 


W°) 2 - 1 

2C 


where a m is the radius of the middle surface of the curved bar given by 


(24) 


(25) 


(26) 


(27) 


Ka-HO -]©“ 

) 


®m — 2 (® "b b) 


(29) 
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For the isotropic case ( k — 1), the equation may be expanded for k — *1 in the following form for h = 1 
(see app.): 


Mr)] 


4 M ( a 2 b 2 b 2 r 2 a' 

fc— i = r — ln- + b 2 ln- + a 2 \n- 
Q \ r l a b r : 


(30) 


where 

Q = (b 2 — a 2 ) 2 — 4a 2 i> 2 (in ^) (31) 

Equation (30) is exactly the same form as the stress equation given by Timoshenko and Goodier (1970) 
for isotropic materials. 


Location of Maximum a' r 

Differentiating equation (26) with respect to r, and setting the resulting derivative to zero, 


SfWM 1 - 0 


the radial location of the maximum <r' may be calculated as 


( T m\ 2k _ ( k+ l )( b / a ) k+1 [(V®)* 1 - 1 
W " (k — 1) [(6/a) i+1 — 1] 

For the isotropic case k— *1, equation (33) is reduced to 



( r _k\ 

b I 2 In b/a 


UJ 

k-n a\(b/a) 2 -l 

where the relationship 




f [( Va )*- 1 - 

l] 1 ln(6/a) fc_1 


! 1 

1 “ Jb - 1 

/ k—+i 


was used. 

The peak radial stress (<r') max when r = r' m may be written as 


(32) 


(33) 


(34) 


(35) 



(36) 


where r' m is given in equation (33). 
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Thin Curved Bar 


When the thickness of the bar is smaU (that is, b— »-a), the expression C given in equation (27) may be 
expanded in terms of t/a (see eq. (16)) as 


(C) b ^ a 



(37) 


in which all the lower order terms up to 0(t/a) 3 cancelled out. 

The expression for given in equation (33) could be expanded in a similar way to yield 



(38) 


Equation (38) shows that the location of (a f r ) max is always on the inward side of the middle surface of the 
curved bar, and that the location of (a' ) max drifts away from the middle surface more slowly than the case 
°f (^r)max ( see eq. (18)) as the value b/a increases. 

The stress equation (36) may be expanded in terms of t/a as 


ha m (b - a) , ] _ (l/2)k(k 2 - l)(t/a) 4 + - ■ • 

. M ^ Jmax |6-* 0 ~ — 2C 


(39) 


Notice that all the terms up to 0(t/a) 3 disappeared. With the application of equation (37), equation (39) 
becomes 


ha m (b — a) 
M 


K), 


Ji.- 


3 

2 


(40) 


which is independent of material anisotropy (that is, the value of k). Notice that the numerical values of 
both equations (20) and (40) are identical. 


DELAMINATION STRESS 


The open-mode maximum delamination stress ctd induced in the curved bar under the end forces P with 
the loading axis offset e (see fig. 1) will be the sum of (<r r ) max due to P without loading axis offset (eq. (15)) 
and (a') max due to M = Pe (eq. (36), see fig. 2): 


or 


&D — (^r)max 4“ (o r r )max 


(41) 




P 

h(b - a) 




(42) 
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where 



For thin-walled curved bar, the values of F\ and F 2 are quite close, or 

Fi « F 2 (45) 

Therefore, the stress contribution from the end moments Pe for the thin-walled curved bar is almost 
proportional to the value of e/a m . 

NUMERICAL RESULTS 

Figure 5 shows the dimensionless delamination stress [6(6 — a)/P](°v)max induced by the end forces P 
plotted as a function of 6/a for different values of anisotropic parameter j3. As b — ► a , the effect of anisotropy 
disappeared, and all the curves converge into a single point giving {[6(6 — a)/I > ](o r r)nuix}& — >. a = 1.5, 
which has been established in equation (20). For low anisotropy 2 < /3 < 5, the delamination stress 
increases monotonically with the increase of 6/a. However, for high anisotropy of Q > 5, the delamination 
stress curves show a slight valley (or dent) in the regions of moderate values of 6/a. In these regions, 
[6(6 - c)/P](< 7 r ) max yields values less than 1.5. Notice that as the value of /? increases, the intensity of 
delamination stress decreases, and that the higher the value of 6/a (that is, the thicker the wall), the higher 
the magnitude of the delamination stress. 

Figure 6 shows the dimensionless radial distance [(r m /a) — 1]/[(6/ a) — 1] of (oy) max point measured from 
the inner boundary of the curved bar plotted as a function of 6/a. It is seen that the effect of anisotropy 
is relatively small and is negligible in the region 6/a < 1.4. As the value of 6/a increases, the location 
of (o’r)max moves away from the middle surface and toward the inner boundary of the curved bar (see 
eq. (18)). 

Figure 7 shows the plots of the dimensionless delamination stress [6a m (6 - a)/M](cr') max induced by 
the end moments M as a function of 6/a. Similar to the previous case, as 6/a— >1, all the stress curves 
converge into one point giving [6a m (6-a)/M](<x') max = 1.5 which was established by equation (40). Notice 
that as the value of k increases, the magnitude of the delamination stress decreases. Figure 8 shows the 
dimensionless radial distance [(r[„/a) - l]/[(6/a) - 1] of (a' ) max point measured from the inner boundary 
of the curved bar plotted as a function of 6/a. The effect of anisotropy turned out to be very small and 
could be neglected in the region 1 < 6/a < 1.7. As 6/a increases (that is, as the wall of the curved bar 
becomes thicker), the location of (cr') max moved inwardly away from the middle surface with a slower rate 
as compared with figure 6 (compare eqs. (18) and (38)). 
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EXAMPLE 

Let the curved bar be made of 25 composite plies with the stacking sequence of [0io/90/0/90/0/90/0io], 
and with the following ply properties: 


E h = 25 x 10 6 lb/in 2 (46) 

E T = 1.2 x 10 6 lb/in 2 (47) 

Gut = 0.75 x 10 6 lb/in 2 (48) 

J'lt = 0.33 (49) 

jp 

= 14T7T = 0.01584 (50) 


where 

Ei, is modulus of elasticity of single ply in fiber direction, 

Ei modulus of elasticity of single ply in the direction transverse to fiber direction, 

Gut shear modulus of single ply, and 
k'LT Poisson’s ratios of single ply. 

Using the mixture rule, the elastic properties of the curved bar composite system may be calculated as 


22 3 

Eg = —.El + —Et = 22.144 x 10 6 lb/in 2 

Zo Z O 

E r = E't - 1.2 x 10 6 lb/in 2 

Gff T = Gut = 0.75 x 10 6 lb/in 2 

22 3 

+ rrt'TL = 0.2923 


(51) 

(52) 

(53) 

(54) 


where 

Eg is modulus of elasticity of laminated composite in 0 direction, 
E r modulus of elasticity of laminated composite in r direction, 
Gg r shear modulus of laminated composite, and 
vg r Poisson’s ratio of laminated composite. 


The two anisotropic parameters (3 and k of the composite system may be calculated as (see eqs (5) 
and (24)) 


/? = 


'1 + |r(l -2vg r )+p- = 6.1799 

tj T C T Q r 


(55) 


k = 


^ = 4.2957 

t, T 


( 56 ) 


It is seen that the composite system under consideration is highly anisotropic. 
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Let the mean radius a m of the curved bar be chosen to be 1 in (ideal specimen size for common fatigue 
test machine, see fig. 1), then the inner and outer radii of the curved bar will be approximately (considering 
ply thickness « 0.0118 in) 


a = 0.085 in ( 57 ) 

b = 1.15 in (58) 

which give 

b/a = 1.3529 in ( 59 ) 

If the loading axis of P has an offset of 

e = 0.375 in ( 60 ) 

( a 1 (VI™*, r m , and r' may be calculated respectively from equations (15), (36), (13), and (33) using 

the numerical values given above. The results are given in table 1, which includes the isotropic case tor 

comparison. 


TABLE 1— PEAK RADIAL STRESSES AND THEIR LOCATIONS 


Items 

Anisotropic 

(/? = 6.1799, k = 4.2957) 

Isotropic 
(/? = fc = i) 

Difference, 

percent 

[h(b - a)/P](cr r ) max 

1.4939 

1.5121 

1.20 

[ha m (b — o)/ P > e](cr r ) max 

1.5093 

1.5198 

0.69 

r m (in) 

0.9767 

0.9774 

0.07 

On) 

0.9814 

0.9812 

0.02 


Notice that the difference in the locations of (<7 r ) max an d (<r') max is 

r' m -r m = 0.9814 - 0.9767 = 0.0047 in 

which is about 40 percent of the ply thickness of 0.0118 in. 

The maximum open-mode delamination stress <td induced by P and M 

<J£) = (<7r)max T (Pr )max 

or 


= — - - f 1.4939 + 1.5093-) 
h(b - a) \ a m J 

(63) 

= „ P x [1.4939 + 0.5660] 
h(b — a) 

(64) 

= - «) 12 05991 

(65) 


(61) 


= Pe may be written as 


( 62 ) 
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The contribution from the end moments M = Pe is about 38 percent, and is roughly proportional to the 
value of e/a m . 


CONCLUDING REMARKS 

Delamination stresses in a semicircular laminated composite curved bar subjected to end forces and end 
moments were calculated, and their radial locations were determined. A family of design curves was 
presented, showing the variation of the intensity of delamination stresses and their radial locations with 
different geometries and different degrees of anisotropy of the curved bar. The information provided can 
be used to select proper geometry of the curved bar fatigue test coupon and to optimize the composite 
stacking sequence so that delamination can initiate at the site of peak delamination stress. 
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APPENDIX— EXPANSION OF STRESS EQUATIONS (26) 

AND (27) FOR ISOTROPIC CASE (k = 1) 

Equations (26) and (27) may be rewritten as (for h — 1) 

( 66 ) 

(67) 

when k — *•!, equations (66) and (67) may be expanded in the neighborhood of k = 1 as follows. 






ifer-i 
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Expansion of o' r 

With equation (72) considered, the stress equation ( 66 ) may be rewritten as 


[ a r( r )]fc — *-i = 


AMa 2 

(k-l)Q 

AMa 2 

(k-l)Q 

AMa 2 
(k - 1)Q 


[«r - ‘] [■ - o'"] - [«r - 

©'-■H-'^reri-ar-e 

er-'lh ©'"]-[© 

-er-en 


r; 

fc-n 


” - ll In Wa)M - 1 


( b/a) 2k — 1 


AMa 2 

(fc-l)Q 




2k 


2k 


- 1 


{bj_ 

(V 


r { b\ k+1 b 

In (* — 1 ) (-) In- 

a \r/ a 

!a) 2k - 1 JJ 


(74) 

(75) 


(76) 


(77) 


The last term of equation (77) may further be expanded as 


K a) 


(fe/g ) fc+1 - l 
( bja) 2k - 1 


l-ener-Hr-c 

-c 


*-l 


£V +1 (*-l)ln£ 

CL J CL 


(78) 

(79) 


Thus equation (77) may be simplified to 


, „ X1 4Ma 2 f[/6\ 2fc 1, r fb\ k+1 fb\ k +\ b\ 

Mr) fc— i = 75 — M- - 1 In - + - - - In — > 

Q ([\a/ J a \rj \aj a J 


or 


(80) 


r / / 4Af 

M r )]*-+1 = -~Q- 


a 2 b 2 , b ,i, r 9 , a 
— 5 - In — + 6 2 In - + a 2 In — 


(81) 


which is the form of the stress equation given by Timoshenko and Goodier (1970) for isotropic curved bar 
under end moments M. 
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Figure 4 • Bending of semicircular 
curved bar by moments at its ends 
(case of figure 2, C). 
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Plots of delamination stress (cr r ) max as a function of b/a for different values of 


Middle surface of curved bar 





Middle surface of curved bar 
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